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ABSTRACT 

Calculations  have  been  performed  to  determine  the  upward 
going  muon  flux  leaving  the  earth's  surface  after  having 
been  produced  by  cosmic  ray  neutrinos  in  the  crust.  Only 
neutrinos  produced  in  the  earth's  atmosphere  are  considered. 
Rates  of  the  order  of  one  per  100  square  meters  per  day  might 
be  expected  if  an  intermediate  boson  exists  and  has  a  mass 
less  them  2  BeV. 
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The  possibility  of  detecting  high  energy  neutrinos  in 

1-3 

cosmic  rays  has  been  discussed  by  many  authors.  In 

2 

particular,  it  has  been  pointed  out  that  an  effective  method 
is  to  use  the  earth's  crust  as  a  target  and  observe  high 
energy  muons  which  are  produced  by  cosmic  ray  neutrinos  coming 
from  the  opposite  side  of  the  earth.  Such  an  experiment, 
if  feasible,  may  well  be  the  only  means  to  yield  direct 
information  on  very  high  energy  neutrino  reactions .  Because 
of  the  rapid  decrease  of  cosmic  ray  neutrino  flux  with 
increasing  energy,  the  number  of  muons  produced  in  this  way 
depends  sensitively  on  the  high  energy  behavior  of  neutrino 
reactions.  In  this  letter,  the  expected  upward  muon  flux 
at  the  surface  of  earth  is  calculated  in  some  detail  under 
certain  assumptions  on  the  nature  of  the  weak  interactions . 

We  first  consider  the  case  that  weak  interactions  are 
transmitted  by  intermediate  bosons,  called  W*,  of  a  mass 
rr^  which  is  not  much  bigger  than  1  BeV.  Since  we  are  only 
interested  in  very  high  energy  neutrinos  (from  about  10  BeV 
to  several  hundred  BeV)  the  dominant  reactions  are  the 
coherent  processes 


and 


V 


+  Z  -  [1~  +  w+  +  z 

(1) 

+  Z  -  p.+  +  w~  +  z 

(2) 

where  Z  is  the  target  nucleus .  Such  a  reaction  may  take 
place  at  a  vertical  depth,  say,  (h  cos  0)  from  the  surface 
of  the  earth  and  produce  a  muon  of  energy  p,  moving  at  an  angle  0 
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with  respect  to  the  upward  vertical  direction.  (See  Fig.  1.) 

Let  u>  be  the  final  muon  energy  after  it  travels  through  the 

distance  h.  At  the  surface  of  the  earth  the  muon  flux  within 

the  energy  interval  du;  and  solid  angle  dQ  is  given  by 

00  ® 

dM  =  doudfi  /  dp.  J  dV  Ny  ( dc/dp.)  (  -dp./ dh)  ~ 1  n  (3) 
w  M-+nV 


where  (da/djj.)  is  the  differential  cross  section  per  atom 
for  reactions  (1)  or  ( 2) ,  n  is  the  number  of  atoms  per  unit 
volume,  V  is  the  energy  of  the  neutrino,  Ny  is  the  neutrino 
flux  per  unit  energy  interval  per  unit  solid  angle,  and 
(~d|j,/dh)  is  the  ionization  loss  of  muon  energy  per  unit  length. 
For  very  high  energy  muons,  the  energy  loss  is  given 

4 

approximately  by  (ft  =  c  *  1) 


-(d^/dh) 


27T  n  Z  aa  (me  1 ) 


(4) 


where  Z  is  the  average  number  of  electrons  per  atom,  a  is 

the  fine  structure  constant  and  m  ,  m  are  respectively  the 

e  n. 

mass  of  electron  and  muon.  For  all  subsequent  calculation 

we  use  silicon  as  the  relevant  atom,  since  almost  all  the 

neutrino  reactions  are  taking  place  in  the  earth's  crust. 

The  neutrino  flux  that  is  produced  in  the  atmosphere 

1  2 

can  be  decomposed  into  three  terms:  ' 


Nv  =  N  (pi)  +  N  ( 77")  +  Nv  (K) 


(5) 


due  to  the  decays  of  p.,  v  and  K,  respectively. 

5 

By  using  the  observed  muon  spectrum  in  cosmic  rays,  it 
is  straight  forward  to  derive  the  total  number 
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D(E^,  cos  9)  dE^  dfi  of  ^.-decays  for  muons  of  energy  E^  which 
decay  in  the  atmosphere,  per  unit  area  per  unit  time.  The 
corresponding  neutrino  flux  is  given  by 

Nv(|i)  S  j  J  D(E^,  cos  0 )  |~  5  -  9(V/E^)2  +  4(v/E^)3"]  (dE^/E^) 

V  L  J 

(6) 

In  (6)  we  include  only  the  mu-neutrinos  from  ^.-decays ,  since 
e-neutrinos  can  not  directly  produce  p,*. 

Using  the  results  of  Greisen^,  the  neutrino  flux  due  to 
TT-decays  is 

00 

V*)  -  3  /  ,(V  [b~'T  cos  e](dW  (7) 

2  V 

3 

where  f(E)  =0.16  (E/BeV)  per  BeV-cm2-sterad  (8) 

and  B_  =  120  BeV. 
t r 

To  obtain  the  neutrino  flux  due  to  K-decays,  we  assume 
the  K-meson  spectrum  has  approximately  the  same  energy 
dependence  as  the  ir-meson  spectrum.  Therefore,  Ny(K)  is 
given  by 

00 

VK)  =  <b  /  "V  [  Bk  T  e“  cos  9  ]  <dVEK>  <9> 

where  BR  =  900  BeV,  f  is  given  by  (8) ,  r  is  the  average 
production  ratio  of  (K/tt)  at  the  same  energy  in  a  high  energy 
nucleon-nucleon  collision,  and  b  is  the  branching  ratio  of 
K^2  decay.  For  definiteness,  we  take  (br)  =  10%. 

Because  of  the  much  greater  fractional  energy  given  to 
the  neutrino  in  a  K^'decay  as  compared  to  that  in  a  7r^2-decay, 
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the  effect  of  Ny(K)  is  especially  important  for  high  energy 
neutrinos.  In  contrast,  the  observed  muon  spectrum  is 
hardly  influenced  at  all  by  the  presence  of  K  2”decay,  Pro_ 
vided  (b  r)  is  small. 

The  complete  expressions  of  the  differential  cross 

sections  for  reactions  (1)  and  (2)  have  been  given  in  the 
6 

literature.  Let  F(q)  be  the  charge  form  factor  of  the 
target  nucleus: 


where  q  is  the  4-momentum  transfer  of  the  nucleus  and  R  is 
its  charge  radius.  For  the  detection  of  cosmic  ray  neutrinos, 
one  is  only  interested  in  relativistic  p  and  W;  further¬ 
more,  (q/p)  and  (qAO  are  both  much  smaller  than  unity. 

(p  and  W  refer  to  the  energy  of  the  muon  and  the  intermediate 
boson  in  the  laboratory  system.)  It  can  be  shown  that  in 
such  a  case  the  differential  cross  section  for  reaction  (1) 

or  ( 2)  with  the  production  of  a  W  of  helicity  s  is  given  by 

3 

do  =  Z2a2  (tt/5)  -1  G(dp/v)  E  c“  Am  (11) 

s  m=0  s  m 

09 

where  Am(x)  =  J  y”24™  B(y)  dY'  (12) 

x 


B(y)  = 


CO 

r 


j 


O  “4  J  ~  o 

(1  +  zz)  (zz  -  yz)  z  dz 


[y"W>]  - 


2»  -1  ,  1 


(l+4y*)  in  y  (1+y  )  -  4  +  (1+y*)  +  £  (1+y  ) 


(13) 
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x  =  (AAO  +  ( 6/p.) 

A  =  2'1  (12) '1/2  R,  (14) 

6  =  2"1  (12)  "1/2  in2  R, 

[L 

G  is  the  Fermi  coupling  constant  and  the  subscripts  s  =  R,  L,  l 

denote,  respectively,  a  right-handed  W±,  a  left-handed  W±, 

and  a  longitudinal  W*.  In  Table  1,  we  tabulate  the  coefficients 

Cm  where 
s 

dx  =  -2  [(2A/W)  +  (64l)J  (1  -  *)  £l  +  (p/vjj  ,  (15) 

d2  =  4  £(2A/W)  +  ( 6/p.)J  2 

and  k  is  related  to  the  magnetic  moment  of  W  by 

(magnetic  moment)  =  ( 2m^)  e  (1  +  <c)  (16) 


It  is  important  to  notice  that  (11)  is  applicable  for 
arbitrary  value  of  (rr^  R/v)  provided  that  (q/|i)  ,  (q AO  are 
small  compared  to  unity  and  that  ji*  and  VT^  are  relativistic. 

Combining  (3) -(11)  the  differential  muon  flux 

2 

(d  M/doidQ)  and  its  integrated  spectrum  (dM/diu)  and  M  at  the 
surface  of  earth  have  been  calculated  numerically  by  assigning 
different  values  of  and  k,  for  cos  9  varying  from  0  to  1 
and  for  w  ranging  from  5  BeV  to  1000  BeV.  Some  typical 
examples  of  (d  M/duudQ)  and  (dM/duu)  are  given  in  Fig.  3  and 
Fig.  2.  The  integrated  values  of  muon  flux  for  uu  a  5  BeV 
are  given  in  Table  II.  This  numerical  calculation  was  made 


with  the  IBM  7090  computer  at  Columbia  University. 

It  should  be  emphasized  that  this  calculated  muon  flux 
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is  only  a  lower  limit  since  it  does  not  include  (i)  muons 
produced  in  the  W±-decays;  (ii)  muons  produced  by  the  non¬ 
coherent  W-producing  processes;  (iii)  neutrino  reactions  in 
which  is  not  produced  in  the  final  state;  ( iv)  neutrinos 
not  produced  in  the  atmosphere.  The  last  effect  has  been 
discussed  by  K.  Greisen.^" 

The  effect  of  (i)  is  important  only  if  the  leptonic 
decays  of  W*  turn  out  to  be  the  dominant  modes.  In  such  a 
case  an  increase  of  muon  flux  &  25%  from  the  above  calcula¬ 
tion  may  be  expected. 

The  effect  of  (ii)  is  important  only  for  relatively  low 
energy  neutrinos.  At  very  high  energy,  the  cross  section  for 
non-coherent  W-producing  reaction  is  expected  to  be  about 
Z  ^  of  that  of  the  coherent  reactions.  The  effect  of  (iii) 
is  due  to  reactions  such  as: 

V  (or  V  )  +  nucleon  -  |x~  (or  p.+)  +  nucleon  +  pions  +  ...  (17) 

r*  M> 

If  W*  does  not  exist  or  if  it  exists  but  its  mass  is 
too  large,  (17)  is,  then,  the  only  relevant  reaction.  Due 
to  the  presence  of  strong  interactions,  the  cross  section 
for  (17)  is  not  known  theoretically.  Nevertheless,  it  is 
expected  to  be  proportional  to  G  (q  )  where  (q  )  is  the 
average  of  the  square  of  the  4-momentum  transfer  given  to 
the  lepton.  Let  (N^)  and  (P^)^  be  the  average  multiplicity 
of  pions  and  the  average  of  the  square  of  the  transverse 
momentum  of  the  pions .  By  assuming  a  random  distribution  of 
the  direction  of  the  vector  (p  )w  we  have 

f~L  TT 

"tot  ~  S2  <q2>  ~  G2  <  V  <pJ!>t 


(18) 
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We  further  assume  that  at  very  high  energy  the  behaviors 
o 

of  (N^)  and  (Pj_>T  are  similar  to  that  in  a  high  energy 
nucleon-nucleon  collision;  i.e., 


and 


constant 


(N  >  CZE® 
t t  c.m. 


where  E  is  total  energy  in  the  center  of  mass  system  and 

c  •  rn  • 

3  is  very  roughly  w  In  such  a  case,  the  total  cross 

2  0 

section  becomes  proportional  to  G  El,  ,  instead  of  the  usual 

c  •  m  • 

2  2 

G  Ec  m  expression  for  a  neutrino  interacting  with  particles 
with  no  strong  interactions .  A  convenient  phenomenological 
formula  for  the  total  cross  section  of  (17)  may,  therefore, 
be  written  as  ( averaged  over  both  neutron  and  proton) 


^ot s  (c/,r)  °2  mp  <v/ys/2 


(19) 


where  m^  is  the  nucleon  mass,  v  is  the  neutrino  energy  in  the 
laboratory  system  (V  »  mp) ,  and  C  is  a  dimensionless  pro¬ 
portionality  constant.  (For  a  high  energy  V  interacting 
with  a  "bare"  nucleon  with  no  strong  interactions  C  =  1  and 
3=2.)  The  resulting  muon  flux  can  be  computed  by  using 
(4)  and  ( 5)  . 

By  carrying  out  these  numerical  calculations,  we  fttid  the 
total  muon  flux  for  uu  s  5  BeV  and  integrated  over  2tt  solid 
angle 

M  =  C  [0.038]  per  (10  meter) ^  -  day  for  3  =  £ 

and 

2 

M  =  C  [0.095]  per  (10  meter)  -  day  for  3=1. 

The  spectra  for  these  two  cases  are  given  in  Fig.  4. 
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Therefore,  it  seems  that  (17)  becomes  important  only  if 

W  does  not  exist,  or  if  it  exists  but  its  mass  is  larger  than 

2 

2  BeV.  For  smaller  values  of  m^,  the  average  (q  )  can  be 
substantially  reduced  by  the 


.  .  2,  21  -1 

1  +  q  /n^ 


factor  in  the  W±  propagator.  Thus  the  muon  flux  due  to 
reaction  (17)  may  be  much  smaller  than  the  above  estimation. 
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TABLE  I 

Table  of  the  coefficients  Cm  (s  =  R,  L  and  l) 

s 

See  (14),  (15),  (16)  for  the  definitions  of  x,  k  ,  d^,  d2 
m  (An.)'1  W2  C™  (AV2)"1  W2|j.  C™  v'^  C 

0  0  0  l(l-K)2  [1+  (n/v)2] 

5-2  +  <c2  4  (1-  K)2(n/v)2  -XV  _1  W  C°  +  dx 

2  I  4(<  -3)x  -4  [2  +  (1  -  /c)  (n/v)  3x  -xd^  +  d2 


3 


a  ** 
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TABLE  II 

- 1  -2 

Muon  flux  in  units  of  (day)  (10  Meter)  ,  integrated 
over  ( 2ir)  solid  angle  and  muon  energy  from  5  BeV  up. 
Only  the  muons  produced  by  reactions  (1)  and  (2)  are 
included. 


\  fC 

(n^/BeV)  \. 

2 

I 

j 

1 

0 

0.5 

6.46 

3.15 

1.48 

0.75 

3.31 

1.67 

0.801 

1 

1.93 

0.998 

0.491 

1.5 

0.833 

0.446 

0.228 

2 

0.432 

0.237 

0.125 
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Fig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


FIGURE  CAPTIONS 


Schematic  drawing  of  muons  produced  by  cosmic  ray 
neutrinos  interacting  with  the  earth's  crust. 


Flux  of  (x-mesons  as  a  function  of  energy  for  k  -  0 
and  different  values. 

Angular  distributions  of  muons  for  =  1  BeV, 

k.  =  0  and  uu  =10  BeV  and  200  BeV.  The  angle  0 
M* 

is  between  the  momentum  of  p.-meson  and  the  verticle 
direction. 


-  Energy  spectrum  of  muons  leaving  the  earth's  surface 
which  have  been  produced  as  the  result  of  reaction  (17) 
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